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The effect of noise is studied in one-dimensional maps undergoing transcriticai, tangent, and 
pitchfork bifurcations. The attractors of the noiseiess map become metastable states in the presence 
of noise, fn the weak-noise limit, a symplectic two-dimensional map is associated with the original 
one-dimensional map. The consequences of their noninvertibility on the phase-space structures are 
discussed. Heteroclinic orbits are identified which play a key role in the determination of the escape 
rates from the metastable states. Near bifurcations, the critical slowing down justifies the use of 
a continuous-time approximation replacing maps by flows, which allows the analytic calculation of 
the escape rates. This method provides the universal scaling behavior of the escape rates at the 
bifurcations. 



I. INTRODUCTION 

At the niacroscale, deterministic dynamics often 
emerges out of the random motion of the many sub- 
units interacting in a complex system. In between the 
micro- and macroscales, the time evolution may admit a 
mesoscopic description in terms of stochastic processes in 
which the deterministic dynamics is affected by noise. In 
such systems, the deterministic dynamics rules the mean 
fields or slow modes emerging out of the many rapidly 
fluctuating degrees of freedom. At the mesoscopic scale, 
these latter can still manifest themselves and perturb the 
deterministic time evolution of the slow modes by fast 
additive or multiplicative noises [TJ [5]. Such stochas- 
tic processes have been considered for Brownian motion 
p], mesoscopic electric circuits and superconducting 
junctions 0, mesoscopic lasers and other optoelectronic 
devices [5] , fluctuating hydrodynamics |7] H] , nucleation 
processes [5], reaction rate theory [TD], nonequilibrium 
chemical kinetics [HI [12] , as well as collective behaviors 
such as synchronization, flocking or swarming in biology 
[I1[I1[T1[T51[I71[TS1[T1]. 

The presence of noise may modify qualitatively the 
time evolution of systems. In particular, noise can ac- 
tivate the crossing of otherwise impenetrable barriers. 
Therefore, states which are stable with respect to the de- 
terministic dynamics may become metastable under the 
effect of noise, as it is the case in nucleation processes or 
chemical reactions [HI [IHl [H] • Instead of remaining for- 
ever in a stable state, the trajectory thus escapes from 
the metastable state with a rate, which corresponds to 
the lifetime of the system in the metastable state. The 
phenomenon of noise-induced escape is known in many 
systems including chemical reactions [TU], nucleation pro- 
cesses [9], Josephson tunnel junctions 0, semiconductor 
lasers [20], or systems with chaotic attractors or fractal 
basin boundaries [23 jH (HIMl [ISIIIS]- 

Analytical methods have been developed to evaluate 
the lifetime of metastable states, especially, in the weak- 
noise limit. In this limit, the lifetime is often observed to 
grow exponentially with the inverse of the noise ampli- 



tude, a behavior which is reminiscent of Arrhenius' law in 
chemical kinetics. In microscopic chemical reactions, the 
noise is of thermal origin and its amplitude is controlled 
by the temperature so that the rate of hopping over the 
energy barrier is proportional to a Boltzmann factor at 
the ambient temperature. In macroscopic complex sys- 
tems, the fluctuations come from the individual motions 
of the multiple active entities interacting in the system 
so that the noise amplitude is related to the size of the 
population of entities and has no thermal origin. In such 
complex systems, the analogue of the activation energy is 
given by the action functional ruling the stochastic pro- 
cess in the weak-noise limit. 

If the noise is weak enough, a variational principle can 
be used to deduce a Lagrangian or Hamiltonian dynam- 
ical system from the action functional, as shown by On- 
sager and Machlup [27], FreidHn and Wentzell [23, and 
others [291 [301 |3l| |32| [33]. Remarkably, this dynam- 
ical system is deterministic and symplectic, i.e., area- 
preserving in problems with one degree of freedom. In 
this respect, the method based on this Hamiltonian dy- 
namics is called the symplectic approach. We notice 
that this approach is similar to the semiclassical method 
of quantum mechanics in the limit where the mechani- 
cal action is much larger than Planck's constant. As in 
quantum mechanics, the action functional plays a key 
role in the evaluation of various properties of interest 

[M1I13IM]. 

In the present paper, our purpose is to use the symplec- 
tic approach to calculate analytically the rate of noise- 
induced escape from attractors undergoing transcriticai 
or pitchfork bifurcations [37'. Such bifurcations occur 
in dissipative deterministic systems controlled by tuning 
parameters. 

At a transcriticai bifurcation, two steady state cross 
each other and exchange their stability. On both sides of 
the bifurcation, the unstable state may form the bound- 
ary with respect to a remote attractor. In this regard, 
this is a boundary bifurcation. In the presence of noise, 
this boundary can be crossed and trajectories may escape 
towards the remote attractor. Near the bifurcation, the 
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proximity of the boundary to the stable state can strongly 
enhance the escape rate in a way to be determined. 

At a pitchfork bifurcation, a stable state becomes un- 
stable, leading to the emergence of a pair of stable states. 
Therefore, the deterministic system becomes bistable 
with two attractors beyond the bifurcation. Noise will 
induce random across the unstable state forming the 
boundary between both attractors. Here also, the escape 
rate from each one of both attractors is deeply modified 
close to the bifurcation. 

As we show in the present paper, the symplectic ap- 
proach allows us to evaluate the effects of these bifurca- 
tions on the lifetimes of the metastable states in the weak- 
noise limit. To investigate these effects, we consider one- 
dimensional maps with additive independent Gaussian 
random variables at each time step. In the weak-noise 
limit, we obtain a symplectic area-preserving map associ- 
ated with the action functional of this Gaussian stochas- 
tic process. Near the bifurcation, we show that an appro- 
priate rescaling of time transforms the two-dimensional 
symplectic map into a continuous-time Hamiltonian sys- 
tem suitable for analytical calculations. In this way, the 
noise-induced escape rate can be evaluated analytically. 

The paper is organized as follow. In Section [IT] we 
briefly present the escape problem for a general one- 
dimensional noisy map. In Section [Till we detail the path 
integral formalism for such noisy systems. Thereafter, 
we perform semiclassical calculation using the steepest- 
descent method in order to obtain the general expression 
of the aforementioned symplectic map. Its phase por- 
traits are analyzed for a specific example. With the help 
of this analysis, we introduce an analytical model in Sec- 
tion |IV| which allows us to calculate an approximate but 
analytic expression for the analogue of the activation en- 
ergy in our context. In Section |V] we present the results 
given by this model for several noisy maps and compare 
them with the results of Monte-Carlo simulations. Fi- 
nally, we give our conclusions in the last section. 



II. NOISY MAPS 

Consider a general one-dimensional map 

Xn+l = f{x„) 



0.5 



(1) 



from the set M of real numbers onto itself. Such maps are 
models of the Poincare recurrence in strongly dissipative 
systems. We notice that the dynamics is deterministic in 
the sense that each initial condition generates a unique 
trajectory: Xn = /"(xq). Typically, the trajectories of 
such deterministic dynamical systems tend to an attrac- 
tor or escape to infinity (see Fig. [T]). The attractors 
of one-dimensional maps can be stationary, periodic, or 
chaotic. The basin of attraction B{T) C M of an attrac- 
tor r is the set of all the initial conditions of trajectories 
converging to F. 

In some circumstances, the deterministic dynamics is 
affected by additive noise and the quantity x is ruled by 
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FIG. 1; Phase space R versus time with typical trajectories of 
the logistic map x„+i = lJ.x„{l — x^) for the parameter value 
jj, = 1.9. The attracting fixed point P of coordinate xr = 
1 ~ and the repulsive fixed point /3 {xp = 0) are plotted 
as straight lines. /3 is one of the two boundaries of the basin of 
attraction -B(r) = [0, 1] of the attractor P. The trajectories 
outside of [0, 1] tend to —oo. Initial conditions xq > \ are 
sent to the region ] — oo, 0[ at the first step. Thereafter, these 
trajectories stay there and tend to — oo. 
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FIG. 2: Same plot as Fig. [Tjbut for the noisy logistic map ([2| 
with = 1.9. For the sake of clarity, each trajectory is plotted 
each thirty steps. Each trajectory starts on the attracting 
fixed point P and fiuctuates for a while before crossing the 
boundary /3 of the basin of attraction B{T) — [0,1]. After 
this crossing, the trajectory is driven to — oo. 



the following noisy map 



(2) 



where {£,n}n=-oo ^ sequence of independent Gaussian 
random variables 





eS„ 



(3) 
(4) 



£ being the amplitude of the noise and (5„„i the Kronecker 
symbol. The trajectories of this system are no longer de- 
terministic since the trajectories {xi}f^Q depends on the 
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sequence of random variables {CiliLo besides the initial 
condition xq. 

One of the consequences is that the noise perturbs the 
stability of the attractors. The boundaries of the basin 
of attraction B{T) can be crossed and the trajectory may 
escape from attractors. Therefore, the attractors T of the 
deterministic dynamics typically become metastable in 
the presence of noise. Noise-induced escape is illustrated 
in Fig. |2]for a noisy map with an otherwise stable steady 
state. 



A. Perron-Frobenius equation and operator 

Although the trajectories of noisy maps are random, a 
deterministic description is recovered for the probability 
density Pn{x) that the system is found in the position x 
at the n^^ time step. Indeed, the probability density is 
ruled by the Perron-Frobenius equation ^56^ given by 

Pn+i{y) = J K{x,y) pn{x)dx (5) 

in terms of a stochastic kernel K(x,y) — g[y — f{x)] 
where g{£,) is the probability density of the noise [5H] . 
For the Gaussian noise ([3|-(|4]) of noise amplitude e, this 
probability density is given by 

In the zero- noise limit e ^ 0, the Gaussian probabil- 
ity density tends to the Dirac delta distribution and the 
Perron-Frobenius equation of the deterministic map is 
recovered [5H] : 

Pn+i{y) ^ J 5[y- f{x)] pn{x) dx (7) 

We notice that, with or without noise, the Perron- 
Frobenius equation is defined by a linear operator 
P, which preserves the non-negativity of the proba- 
bility density as well as its normalization to unity: 
J Pn{x) dx = 1. In order to analyze the long-time be- 
havior of the probability density, an eigenvalue problem 
can be posed for the Perron-Frobenius operator 

-P*a(a;) = Xa*a(a;) (8) 

Since the probability density is not expected to grow 
without bound, the eigenvalues should belong to the unit 
disk: IxqI < 1. 

The border of the unit disk contains the single eigen- 
value xq = 1 a the probability density converges to a 
unique stationary solution of the Perron-Frobenius equa- 
tion. In this case, all the other eigenvalues are strictly 
inside the unit disk, |xa| < 1 for a ^ 0. As a conse- 
quence, the components of the probability density cor- 
responding to these other eigenvalues are exponentially 
damped in the long-time limit, allowing the convergence 



of the probability density to the eigenfunction "^q asso- 
ciated with the leading eigenvalue xo = 1- 

In contrast, the probability density is expected to van- 
ish if trajectories escape to infinity as observed in Fig. |2] 
In this further case, the border of the unit disk does not 
contain any eigenvalue, otherwise the asymptotic proba- 
bility density would not be vanishing. Instead, the long- 
time behavior is dominated by the eigenvalue which is 
the closest to the border of the unit disk: xo < 1 such 
that Ixal < Ixol for a ^ 0. After n iterates of the Perron- 
Frobenius operator, the probability density is damped by 
a factor Xo so that we may expect an exponential decay 
in this situation. The corresponding eigenfunction ^'o(a:) 
is non-negative according to the Perron-Frobenius theo- 
rem. 



B. Escape rate and mean escape time 

Thanks to the probability density Pnix), we can write 
down the probability that the escape from the attractor 
r occurs at a time greater than n in the noisy system Q : 

P{{t>n}) = Gn= f pn{x)dx (9) 
JBir) 

where B{T) is the basin of attraction of F. This proba- 
bility is also called the survival probability. On the other 
hand, the probability to escape at a time n is given by 

P({< - n}) = P„ = G„-i - G„ (10) 

We notice that the normalization condition X^^i P"- ~ 
Go = 1 is satisfied. 

If trajectories escape to infinity, the escape rate 7 of 
the attractor of F is defined as the long-time limit of the 
exponential decay rate of the survival probability: 

7= lim --lnG„ (11) 

n — »oo Ti 

In this case, the escape can alternatively be defined 
in terms of the leading eigenvalue xo of the Perron- 
Frobenius operator. Since the probability density is 
damped as Xo for n — > 00, so does the survival prob- 
ability ([9]) and we infer that 

7 = -lnxo (12) 

The inverse of the escape rate is the lifetime of the 
metastable state corresponding to the eigenfunction 

r^^ (13) 

7 

On the other hand, the mean escape time can be de- 
fined as the statistical average of the random time the 
trajectory escapes from the basin -B(r): 

00 00 
(t)^^nP„-^G„ (14) 

n— 1 n— 
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If the decay is exponential in the long-time limit, G„ ~ 
e""*"", the mean escape time is given by (t) ~ (1 — e"''')"^. 
In the weak-noise limit where the escape rate is expected 
to vanish as 



7 ~ exp 



Wo 

e 



(15) 



for e ^ 0, the mean escape time can be related to the 
escape rate itself by the simple relation 



lim 7 (i) = 1 

7— >o 



(16) 



Otherwise, the mean escape time is not simply given as 
the inverse of the escape rate because, in general, the 
mean escape time depends on the transients before the 
long-time exponential decay at the escape rate. 



the present process by Eq. (14 1. If the attractor F of 



the deterministic dynamics is strictly inside its basin of 
attraction B{T), we expect that, in the weak-noise limit, 
the relaxation rate should behave as 7 ~ exp{—Wo/e), 
for e ^ 0, and be related to the mean escape time (t) 
and the corresponding escape rate 7 according to 



lim7(i> 

7— 



lim7(t) 
7— ►o 



1 



(19) 



For both escape and relaxation, the determination of 
the lifetime of the metastable state thus goes by the eval- 
uation of the quantity Wq which is the analogue of the 
activation energy in systems with thermal noise. 



III. PATH INTEGRALS AND THE 
SYMPLECTIC APPROACH 



C. Relaxation rate and mean first-exit time 

If trajectories do not escape to infinity but randomly 
jump between two attractors of the deterministic dynam- 
ics, the probability density is expected to converge to- 
wards a stationary density given by the eigenfunction 
^'o(a:) of the Perron- Frobenius operator corresponding 
to the eigenvalue xo = 1- In this case, the escape rate 
(111 vanishes and the quantity of interest is given by the 
next-to- leading eigenvalue xi such that |xa| < Xi < 1 for 
a ^ 0, 1 if this latter is unique. The eigenvalue xi con- 
trols the convergence of the probability density towards 
the stationary state. This convergence is exponential, 
which defines the relaxation rate 



7 = -Inxi 



(17) 



This relaxation rate can be interpreted as the rate of 
jumps between the two attractors due to the noise. In 
some cases, this relaxation rate can be defined in the 
presence of escape if there is a clear separation of time 
scales between both. 

A mean first-exit time can also be defined as the sta- 
tistical average of the random time of first exit from a 
given basin of attraction B(T). In this scheme, the tra- 
jectories starts with an initial condition inside the basin 
B(T), evolves in time according to the noisy map (U, and 
stops as soon as the trajectory escapes from the basin 
B{T). The Perron-Frobenius operator Q of this first-exit 
process is defined by restricting the domain of integration 
to the basin B{T): 



Pn+iiy) = 



B(T) 



K{x,y) pnix) dx 



(18) 



For a bounded basin B{T), all the eigenvalues of this 
further operator are strictly inside the unit disk and an 
escape rate 7 can be defined by Eq. ( 12 ) in terms of the 



leading eigenvalue xo ; in the previous Subsection |IIB[ 
Furthermore, a mean escape time (t) can be defined for 



The aim of the present section is to obtain the sym- 
plectic map ruling the process in the weak-noise limit. 



A. Path integrals 

Starting with an initial density po and iterating Eq. 
([5]), one can obtain the probability density at time n in 
the form 



Pn{Xn)= dxoKn{xo,Xn) Poixo) 



(20) 



where Kn is an integral kernel resulting from n iterations. 
This kernel is the propagator associated with the time 
evolution of the stochastic process. For the Gaussian 
noise ([6]), this kernel is given by the multiple or path 
integral 



Kn{xa,Xn) = 



1 



dx,, 



da;i 



(21) 



X exp 



1 



Wn{xo . . . ^ Xn—±^ Xyi) 



where 



n-l 



Wnixo, 



1 " " 

") = 9 E [^'+1 - /(^')]' (22) 



i=0 



which defines the action functional of the Gaussian 
stochastic process. 



B. Symplectic map 



In the limit of weak noise e ^ 0, the propagator 



(21 1 can be integrated by the steepest-descent method. 



which will select the extrema of the action functional 



( 22 1 . These extrema are obtained by applying Hamil- 



ton's variational principle for the action function (22 1 
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The selected paths are thus the trajectories of a sym- 
plectic deterministic system [39] . With this Hamiltonian 
formalism in mind, the functional (22 1 can be written in 



terms of a generating function F of the first kind 



W„{xo, ...,Xn) 



ji-i 

E 

1=0 



F{xi,Xi+i) 



where 



1 2 

F{x^,X^+l) = - [Xi+i - f{x^)] 



(23) 



(24) 



Introducing the momentum pi canonically conjugated to 
the position Xi, the evolution over one time step is con- 
ceived as a canonical transformation corresponding to 
this generating function 



Pi+i 



Pi 



dF(xi,Xij^i) 
dxi 



(25) 



With the generating function (24), we obtain the two- 
dimensional map 



,+1 = !{xi) + j^^ 

— Pi 



(26) 



This map is symplectic because it preserves the symplec- 
tic differential two-form 



dxi+i A dpi+i = dxi A dpi 



(27) 



where A denotes the exterior or wedge product |41j . Ac- 
cordingly, the Jacobian of the map takes the unit value, 
det D(j) = 1, so that the map is area-preserving in the two- 
dimensional phase space M = {{x,p) : x £ R,p e M}. 

Such area-preserving maps are expected to be typically 
chaotic as it is the case for the standard map [411. The 
multitude of possible trajectories reflects the complexity 
of the landscape of the multidimensional functional ( [23| . 
This landscape unveils many local extrema among which 
the local minima should be selected in order to guaran- 
tee the local normalization of the probability distribution 



underlying the propagator (21) 



We emphasize that the symplectic map ( 26 1 has special 



properties which are specific to the problem at hand. 

First, the map (26) leaves invariant the subspace 
p = of the phase space A4. In this subspace, the 
two-dimensional map ( 26 1 reduces to the noiseless one- 
dimensional dissipative map x^+i = /(x;). In conse- 
quence, the attracting fixed points of the original map 
/ are hyperbolic fixed points of the symplectic map ( |26[ ) 
that lie in the subspace p — 0. Note that the repulsive 
fixed points of the map / are also hyperbolic fixed points 
of (26) and also lie on the subspace p — 0. 



Secondly, the symplectic map ( 26 1 is not uniquely in- 
vertible and may admit several inverses. This is the 



consequence of the noninvertible character of the one- 
dimensional map / itself, which often has extrema where 
its derivative f'{x) vanishes. Indeed, the inverses of the 
two-dimensional map ( 26 1 are directly associated with 



the inverses /-I of the one-dimensional map /, such that 
fi^ [/(a;)] = X for j ~ l,2,...,m. The inverses of the 



two-dimensional map ( 26 1 can thus be written as 



= ^{Xt+l -Pi+i) 



Pi 



Pt 



i+l 



Pi+l)] 



(28) 



with j — 1, 2, m. The inverses of one-dimensional map 
become singular at the positions such that f'{xk) = 
with fc=l,2,...,m. These loci correspond to the borders 



of the domains where the inverses ( 28 1 are defined 



X -p = f{xk) 
with k = 1, 2, ?7i 



such that 



f{xk) = (29) 



For noninvertible two-dimensional 



maps such as Eq. ( 26 ) with a vanishing denominator, the 



concepts of focal point and prefocal curves (or lines) have 
been introduced [42j i43j . They will be defined here below 
in the discussion of the phase-space structures generated 
by these maps, which are different in this regard from 
those of the invertible area-preserving maps such as the 
standard map [41] . 



C. Phase space 



This section is devoted to the phase-space structures 



generated by two-dimensional symplectic maps ( 26 1 as 



sociated with one-dimensional maps / undergoing bifur- 
cations. The study of typical phase portraits of these 
symplectic maps is important because the contributions 
to the propagator (21) are evaluated in the weak-noise 
limit in terms of the trajectories of these maps. More- 
over, statistical averages are related to the propagator 
and, as we shall see, the study of phase space gives an 
insight on their properties. 



1 . Invariant subspace p = 

As mentioned here above, the straight line p — is an 
invariant subspace of the phase space M — {(x,p) : x € 
M.,p e M}. This line is the a;-axis where the dynamics 
is ruled by the one-dimensional map x^+i = f{xi) of the 
noiseless macroscopic system. The trajectories outside 
the invariant subspace are specifically associated with 
noise. This invariant subspace divides the phase space 
in two half planes and, therefore, constitutes a major or- 
ganizing geometric structure. 



2. Fixed points 



The symplectic map (26 1 has two types of fixed points. 
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The fixed points x = f{x) of the one-dimensional map, 
which are found in the invariant subspace p = Q. For 
these fixed points, the linearized map has the eigenval- 
ues A+ = f'{x) and A_ = f'{x)~^. Accordingly, these 
fixed points are hyperbolic (or possibly parabolic in the 
cases of marginal stability). If the point is attracting 
(resp. repelling) for the one-dimensional map, the unsta- 
ble (resp. stable) direction is transverse to the invariant 
subspace p= 0. 

On the other hand, the symplectic map may admit 
genuine fixed points outside the invariant subspace p = 
for 



(30) 



These fixed points have eigenvalues given by the roots of 
the characteristic equation: 




+ [pr ix)-2]A+l = 



(31) 



so that they can be cither elliptic or hyperbolic (or also 
parabolic in the marginal cases). 



3. Singular lines, focal points, and prefocal lines 

Because of its vanishing denominators, the two- 
dimensional map (26 1 is not defined on the singular lines 



5s,k^{{xk,p)- f{xk)^0,p£M.} 



with k = l,2,...,m. The implications of these singu- 
lar lines on the phase-space structures generated by the 
map have been discovered by Mira and coworkers who 
introduced the concepts of focal points and prefocal lines 

nana]. 

As long as the numerator corresponding to the denom- 
inator is not vanishing, the points of the singular lines are 
sent to infinity by the map (26 1. However, this might not 



be the case if the numerator vanishes with the denomi- 
nator at the so-called focal points HH US] 



Zk = ixk,p = 0) 



(33) 



where ^ is the parameter and a the slope of the line at 
the focal point. The image of this line is given by 



X' = /(Xfe+0 + 



Hi) 



0(0 



o{0 



(35) 

where we have expanded in powers of the parameter ^ 
around ^ = 0, using f'{xk) = 0, and assuming that 
f"{xk) =/= 0. The image 0(7) is a line for varying values 
of the parameter ^. Taking the limit ^ = 0, we obtain 
the image of the focal point: 



HZk) 



(36) 



p = 



We notice that the image depends on the slope a given to 
the line 7 in Eq. ( 34 1 , whereupon the images are multiple 



and form a line called the prefocal line: 
■■ x' =p' + f{xk) 



(37) 



which coincides with the border ( |29[ ) of the domains of 
definition of the inverse two-dimensional maps ( 28 1 . Ac- 



cordingly, the prefocal line is mapped onto the focal point 
under the inverse map ( 28 ) having the prefocal line in its 
domain of definition: 



(t>j ^ ((^zj = Zk 



(38) 



In this sense, the prefocal line constitutes the image of 
(32) the focal point, Sz^ — 4'{zk), showing that a point can 



be mapped onto a line in such two-dimensional map ( 26 1 
with a vanishing denominator |42l 143] . 



4-. Global stable set and unstable manifold 

The structures generated by the map <j) in its phase 
space M. can be analyzed in terms of the stable and un- 
stable sets associated with the saddle points. 

The unstable set 'W^{z) of the hyperbolic fixed point 
z — {x,p) consists in the set of points that converge to z 
under backward iterations of 4>. This can be expressed in 
terms of the union of the successive images of the local 
unstable manifold Wj"j,(2;), i. e., the unstable set in a 



Since the map ((26| displays an undetermined ratio 0/0 neighborhood of z: 



at each focal point, the image of the focal point is not 
necessarily sent at infinity and we may wonder what is 
the phase-space structure corresponding to this image. 
To answer this question, let us consider a line 7 going 
through the focal point and defined by the following para- 
metric equations: 



(39) 



X = Xk+ ^ 



(34) 



Even if is not invertible, the images of the local man- 
ifold Wio(.(2) will be uniquely determined. Accordingly, 
the set W"(2) is a manifold in the phase space and is 
called the unstable manifold [44 . 

On the other hand, the stable set of the hyperbolic 
fixed point z is defined as the set of points that converge 
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FIG. 3: The phase space A4 of the symplectic map (261 
associated with the exponential map = f^Xn exp{—Xn} 

with fi = 2.6. Global unstable set (in red) and global sta- 
ble manifold (in black) of respectively the fixed point F and 
the boundary l3. At the intersections of the two sets, the 
heteroclinic trajectory (black cross) between these two fixed 
points. In light blue, some bounded trajectories inside the 
region delimited by the two sets. 



to z under forward iterations of (/>. It can be obtained 
as the union of successive pre-images given by the possi- 
bly multiple inverse maps <j)~^ acting on the local stable 



manifold Wi^^iz) 



(40) 



If (j) is invertible, this set is a manifold and one speaks 
about the stable manifold of z. However, if there are 
multiple inverse maps, the global stable set may consist 
of disjoint pieces and is not a manifold [44, ,45^. Since we 
here consider noninvertible maps /, the symplectic map 
( |26[ ) is noninvertible so that we speak about the global 
stable set (instead of stable manifold). 

At the intersection of the global stable set and the 
unstable manifold associated with one and the same hy- 
perbolic fixed point z, we find the so-called homoclinic 
orbit W%z)nW{z). 

If the global stable set and the unstable manifold are 
attached to different fixed points, we are in the presence 
of a so-called heteroclinic orbit ^^{z) n yV"(2;'), which 
plays an important role in the following. 

As an example, a typical phase portrait of the sym- 
plectic map ( 26 1 is depicted in Fig. [s] for the expo- 
nential map Xn+i — Ma;n exp(— x„). For the noisy ex- 
ponential map with the selected value of the parameter 
the escape scenario is the same as the one detailed 
in Fig. [2] for the logistic map. Starting from the at- 
tractive fixed point x-p — In /x, trajectories escape after 
crossing of the boundary x^j = 0. Once this boundary 



is crossed, the trajectory tends to —00. We notice that 
the fixed points of the one-dimensional map correspond 
to the fixed points F = (In/i, 0) and /3 = (0,0) of the 
two-dimensional map. Moreover, the stable set yV"^(F) 
and the unstable set yV^(/3) are contained in the invari- 
ant subspace p = 0. The construction of the unstable 
manifold yV"(F) and the global stable set yV"'(/3) seen in 
Fig. [3] have first required the determination of the mani- 
folds directly emanating from the hyperbolic fixed points. 
This has been performed with an algorithm specially de- 
voted to noninvertible maps and called the search circle 
algorithm [3U|3B]. Thereafter, we were able to compute 
significant parts of the global invariant sets using back- 
ward iterations with <\)^^ for the global stable set W(/9) 
and forward iterations with (\) for the unstable manifold 
yV"(r). We have also computed with precision the het- 
eroclinic orbit W(/3) n yV"(r) using an algorithm de- 
scribed in Ref. [47 . An attentive examination of Fig. [s] 
shows that the global stable set yV"'(/3) of the boundary 
fixed point /3 forms closed loops around the focal point 
z — (1,0). This feature is inherent to the map (26 1 which 
is not defined on the singular lines where the derivative 
of the one-dimensional map / vanishes. For the expo- 
nential map, there is such a singular line at a: = 1 and its 
intersection with the subspace p = where there is a 0/0 
limit gives the focal point z — (1, 0) where /'(I) — 0. The 
closed loops of the global stable set yV"'(/3) are attached to 
this focal point. This feature is related to the noninvert- 
ibility of / for which the focal points divide the configu- 
ration space X into different pre-image regions. Since the 
nondefinition sets are lines that divide the phase space 
into distinct regions, the only way for a curve that is for- 
ward invariant to cross these lines is to pass through focal 
points, as here observed for the global stable set >V(/?). 
We refer the reader to the papers [451 |33J [IB] for more 
information about this issue. 



5. Dynamics close to the bifurcation 

Close to the bifurcation, the dynamics undergoes a 
critical slowing down. As a consequence, the symplec- 
tic map ( 26 1 performs small steps at each iteration in 



the phase-space domain where the fixed points bifurcate. 
It turns out that the dynamics of the map is similar to 
a continuous-time fiow. In order to illustrate this phe- 
nomenon, we continue to consider the exponential map 
Xn+i = /iXn exp(— a;„) and we depict in Fig. |4]its phase 
portrait close to the transcritical bifurcation it under- 
goes at the parameter value /i = 1. At this bifurcation, 
the two fixed points xp — and x-p = In /i cross each 
other and exchange their stability. The unstable of these 
fixed points is the boundary of the basin of attraction 
of the other so that a boundary bifurcation happens at 
the transcritical bifurcation. Below and above the bi- 
furcation, an elliptic island exists in the lower half plane 
p < around the extra fixed point approximately located 
at X ~ {^1 — l)/2 and p ~ —{^ — 1)^/4. At the bifurca- 
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tion, the three fixed points coalesce at the origin with the 
surrounding island. 

Figures |4]^a) and (f) depicts the phase portraits for 
values of fi far from the bifurcation where we observe that 
the global stable set and the unstable manifolds are well 
developed. The trajectories in light blue are trapped and 
form an elliptic island composed of invariant circles and 
chaotic zones of moderate extension. At these values, the 



another fixed point. This means that the series (41 1 con- 



system ( 26 1 is clearly not close to integrability since the 



global stable set and the unstable manifold form typical 
heteroclinic structures. 

For values of /i closer to the bifurcation, Figs. |4|b) 
and (c) show that most of the trajectories of the elliptic 
islands are pretty regular lines bounded by the global 
stable set and the unstable manifold nearly forming a 
separatrix. The chaotic zones have disappeared on tiny 
scales in phase space. Hence, we can conclude that, close 



to the bifurcation, the system (26 1 tends to complete 
integrability. 

It is interesting to understand how the transition oc- 
curs between the two situations. In Fig. |4j^d), the phase 
space is depicted for an intermediate value between these 
of Figs.|4|c) and (f). Although the phase portrait is sim- 
ilar to the one in Fig. |4|^c). Fig. |4|^e) shows a zoom of 
the global stable set and the unstable manifold around 
the origin revealing heteroclinic structures typical of non- 
integrable dynamics. 

We conclude that, as the bifurcation is approached, 
the slowing down of the dynamics tends to transform 
the discrete-time map into a continuous-time flow with 
one degree of freedom which is thus integrable at the 
bifurcation. This idea will be elaborated upon in the 
next section. 



tains nonlinear terms. Let us suppose that Cfc with k > 1 
is the first non- vanishing coefficient besides Ci. 

A bifurcation occurs if the fixed point x — looses or 
gains stability. This is the case at the value of fi where the 
linear stability eigenvalue of the map / at a: = is equal 
to unity. Since the linear stability of a; = is /'(O) = 
Ci, the bifurcation occurs at the critical parameter value 
/i = /ic where ci(/ic) — 1- Near the bifurcation, the 
dynamics of the map is slowed down since its slope is 
close to unity, ci ~ 1, which limits the iterations to small 
steps. Consistently, the other fixed points of the one- 
dimensional map are approximately given by x ~ [(1 — 
ci) / Ck]'^ ^ , which coalesce onto the origin x = at the 
bifurcation. 

In this regard, we can suppose that trajectories 
{xn}t=~oo of the map can be seen as trajectories x{t) 
of a continuous-time stochastic dynamical system with 
the time t = n. Since the time step is unity. At = 1, 
differences become derivatives: 



Xn+i - Xn= = [e^'a;(<) - x{t)\ ^ 
. 1 .. 



(42) 



The critical slowing down near the bifurcation justifies 
the truncation of this expansion to the term with the 
first derivative with respect to the continuous time t = n, 
which amounts to assimilate the difference between two 
successive time steps with a time derivative. 

To implement this idea, the position x„ at the previ- 
ous step is subtracted from the noisy map ([2|, which is 
rewritten as 



Xn+l- Xn = f{Xn) - Xn+ £.n (43) 



IV. CONTINUOUS-TIME LIMIT AT 

BIFURCATION Introducing the function 



A. Hamiltonian flow 

In this section, we establish the correspondence be- 
tween the symplectic map ( |26[ ) and a continuous-time 
flow with one degree of freedom in the neighborhood of 
the bifurcation. This correspondence provides us with 
a method to compute analytically the properties of the 
noisy map close to the bifurcation and, in particular, the 
rate of noise-induced escape from the attractor. 

With this aim, we consider a general map of the form: 

oo 

}{x) = Y,c,x' (41) 

i=0 

where the coefficients Ci(/Lt) depend on the parameter /i. 
For the transcritical and pitchfork bifurcations, we take 
Co = which implies that a; = remains a fixed point 
across the bifurcation. For a bifurcation to happen, we 
must moreover suppose that the map / possesses at least 



g{x) = fix) - X (44) 

the noisy map can be replaced by the Ito stochastic dif- 
ferential equation 

^=g{x)+r^{t) (45) 
with the Gaussian white noise 

m) = (46) 
ivitUt')) = e5{t-t') (47) 

We notice that the time integral of this noise defines the 
Wiener process W(t) = jlri{t')dt'. Weemphas ize the re- 
placement of the noisy map by the stochastic differential 
equation is justified by the critical slowing down of the 
dynamics close to the bifurcation. 

In the weak-noise limit, this stochastic process can be 
analyzed thanks to the theory by Onsager and Machlup 
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FIG. 4: Phase portraits of the symplectic map ( |26[ ) associated with the exponential map x„+i = fiXn exp{~Xn) for different 
values of the parameter fi close to the transcritical bifurcation a.t ^ — 1: (a) fi = 0.1; (b) /i = 0.5; (c) fi = 1.2; (d) /i — 2.0; (f) 
= 2.4. The panel (e) is a zoom of (d) around the fixed point (0, 0). 
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[27] or Freidlin and Wentzell US] in terms of a Hamilto- 
nian flow with one degree of freedom. Its Hamiltonian 
function is given by 



H 



■ + g{x) p 



(48) 



This Hamiltonian function rules the Hamilton-Jacobi 
equation derived from the Fokker-Planck equation asso- 
ciated to the equation (45 1. The trajectories are now 



solutions of Hamilton's equations 

X = g{x)+p 
P = -9'{x)p 



(49) 



where p is the momentum variable canonically conju- 
gated to the position x. We notice that, as for the sym- 
plectic map (26 1, the fixed points of the flow x = g(x) are 



hyperbolic fixed points of Eqs. (49 1 that lie on the invari- 
ant subspace p = of the phase space. At least, one of 
the fixed point xy of the fiow x = g{x) is an attractor. Its 
basin of attraction may be limited by another fixed point 
xp which is unstable. For the exponential map, the phase 
portraits of this Hamiltonian flow and the corresponding 
symplectic map are depicted in Fig. [5] 



The heteroclinic orbit and its action 



For the stochastic process (45), the attractor xt is 
found in a well of the kinetic potential deflned by [49l [50] 



U{x) = — j g{x)dx 



- / f{x)dx 



(50) 



On the other hand, the unstable flxed point xp is the top 
of a barrier for the noise-induced escape from the poten- 
tial well. The escape over this barrier is an activated 
process with an escape rate given by Eq. (15) where the 



constant Wq is the action of a special trajectory of the 
Hamiltonian flow (49 1. This trajectory is given by the 



separatrix connecting the hyperbolic flxed point F corre- 
sponding to the attractor to the other hyperbolic flxed 
point /3 at the top of the barrier [i^l (SHI [5T]- See Fig. 
|5ja) . The parametric equation of this separatrix is 



p{x) = -2g{x) 



(51) 



and the action of this solution between the hyperbolic 
points (a;r,0) and (a;/3,0) is given by 



p{x) dx 



(52) 



where the superscript 'c' stands to recall that this is the 
action of the continuous-time Hamiltonian system. This 
action yield the activation barrier for the escape problem 
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FIG. 6: Evo lution of the heteroclinic orbit of the symplec- 
tic map (26 1 associated with the exponential map x„+i = 
/ia;„ exp(— x„) as the parameter jj, varies from 0.4 to 2. 
Around the boundary bifurcation, the heteroclinic orbit tends 
to a continuous curve which is the separatrix of the Hamilto- 
nian flow ( 49 1 . 



in the stochastic process ( 45 1 from the attracting flxed 
point xr through the boundary x^, as if the crossing 
by the point xp was the only possible channel of escape 
from xr ■ The integration of the function g gives the flnal 
result: 



= 2 [U{xp) - U{xr)] 



(53) 



in terms of the kinetic potential ( 50 1 



We notice that the separatrix of the Hamiltonian flow 



(49) is the hmit of the heterochnic orbit h = W(/3) n 
yV"(F) between the global stable set of the boundary 
point (3 and the unstable manifold of the attractor F. 
Therefore, we expect that the action (52) or (53) of the 



separatrix is the limit of the action of the heteroclinic 
orbit of the symplectic map evaluated by the sum: 



1 



+CX2 



71 — — oo 



(54) 



The approximation of the map action (54) by the flow 



action ( 53 1 is justifled as the bifurcation is approached 
as seen in Fig. [6] where we observe that the heteroclinic 
orbit tends to the separatrix of the Hamiltonian flow (49 1. 



V. RESULTS 

To investigate the relation between the activation bar- 
rier Wo and the actions ( 53 1 or ( 54 ) for different types 
of bifurcations, we have carried out a Monte-Carlo sim- 
ulation generating 2 x 10^ escape events starting from 
the attractor F. We have computed the distributions of 
escape times for different noise amplitudes e. As these 
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distributions slowly decrease exponentially in time in the 
limit e 0, the best estimator of their average is simply 
the statistical average ^T^ according to Eq. (TE^. In this 



way, we have obtained the mean escape time ( 14 1 as a 
function of e. For s —^ 0, this function shows exponential 
increase according to Eq. (IT5|). By fitting, we have com- 



puted the activation barrier Wq. It is important to notice 
that Eq. (15 1 is valid under the condition e << Wq on 



the values of the noise amplitude we had to consider. 

In the case the one-dimensional map has several at- 
tractors, we have to consider the first-exit process from 
one basin of attraction as described in Subsection III CI 



A. Transcritical bifurcation 

In this section, we study two simple maps which un- 
dergo a transcritical bifurcation. 

1. Logistic map 



We first consider the logistic map 



(55) 



which is a well-known model in the fields of hydrody- 
namics [52j [53] , glass formation dynamics [54] , and pop- 
ulation dynamics [371 [5S]. This map possesses two fixed 
points : 



xi ^ 

X2 = 1 - l//i 



(56) 
(57) 



The logistic map undergoes a transcritical bifurcation at 
fj- — 1 where these two fixed points exchange their stabil- 
ity. For /J, < 1, a:i is stable while X2 is unstable and forms 
with its non-trivial pre-image f~^{x2) the boundaries of 
the basin of attraction of xi: B{xi) =]x2, f~^{x2)[- For 
/X > 1, X2 is stable while is unstable and forms with 
its non-trivial pre-image the boundaries the basin of at- 
traction of X2: B{x2) =]0, 1[. In both cases, the point at 
infinity is a second possible attractor for initial conditions 
outside the described basin of attraction. 

For the noisy logistic map, the escape to infinity may 
occur via the two boundaries of the basin of attrac- 
tion B{T). Both boundaries correspond to the same 
fixed point xp since one of the boundaries is the non- 
trivial pre-image f^^{xp) of the other. The phase por- 
trait of the symplectic map ( 26 1 is depicted in Fig. [7] 
where we observe similar structures as in Fig. [3] for the 
exponential map. In particular, an elliptic island sur- 
rounds the non-trivial fixed point at x = (fi — 1)/ (2/i) 
and p = —(/X — l)^/(4/Lt). The symplectic map (26 1 has 
two inverse maps which are defined in the domain 
p> X — ifi/i) limited by the prefocal line (37). The map 
( 26 1 is not defined on the singular line a; = 1/2 where 
fjx) = and a focal point exists at 2; = (1/2, 0). This fo- 
cal point is the passage of the global stable set yV^(/3), as 




FIG. 7: Phase portrai t of the symplectic map ( 26 1 associated 
with the logistic map ( 55 I for p : 



1.9. 



Under this condition, 
the attractor corresponds to the fixed point F = (2:2,0) with 
Eq. ( |57[ ). The fixed point of the escape barrier is /3 = (0,0). 
Note that manifolds of the global stable set have to pass by 
the focal point z = (1/2,0) to cross the nondefinition set 
given by the singular line {{1/2, p) : p £ R}. The heteroclinic 
orbit h = W°(/3) n W"(r) is depicted by crosses. 




FIG. 8: Action of the heteroclinic orbit Wh and activation 
barrier Wo (with error bars) from the Monte-Carlo simulation 
versus the parameter value ^ for the logistic map (551. The 



blue hatched curve is the analytical action VVJ calculated near 
the transcritical bifurcation and given by Eq. ( 58 1 . 



seen in Fig.[7] The heteroclinic orbit h = W(/3)n>V"(F) 
connecting the fixed points F to /? was constructed with 
an algorithm described in Ref. [47j. This heteroclinic 
orbit effectively controls the process of noise-induced es- 
cape to infinity of trajectories issued from the attractor. 
Monte-Carlo simulations have allowed us to compute the 
escape rate and the activation barrier Wq, verifying in 
Fig. [8] that it coincides with the action of the heteroclinic 
orbit. Near the transcritical bifurcation, the continuous- 
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time approximation provides us with the following ana- 
lytical expression for this action: 



(58) 



obtained from Eq. ( 53 1 with the potential ( 50 1 . This for- 



mula applies to both sides of the transcritical bifurcation 
at /X = 1 because of the absolute value | • |. Indeed, the 



fixed points ( 56 1 and ( 57 1 exchange their role of attractor 
and barrier at the bifurcation. Figure |8] shows the excel- 
lent agreement between the activation barrier Wq and Eq. 
( 58 1 as the transcritical bifurcation is approached in the 
limit 1. At the bifurcation, the activation barrier 

vanishes as the cube of the control parameter = ^—1. 
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2. Exponential map 



The exponential map 

a^n+i = M a^ri exp (— a;„) with /i > (59) 

shares several features with the logistic map. It has two 
fixed points: 



xi = 
X2 — In /X 



(60) 
(61) 



which cross each other in a transcritical bifurcation at 
/i = 1. Moreover, the point at a; = — oo is also attracting. 
For < /i < 1, xi is stable, X2 is unstable, and the basin 
of attraction of xi is B(xi) =]a::2, -l-oo[. For fi > I, X2 is 
stable, xi is unstable, and the basin of attraction of X2 
is B{x2) =]xi,-|-oo[. Contrary to the logistic map, these 
basins of attraction are not compact and extend up to 
X = +00 where the exponential map vanishes and sends 
the points toward the fixed point xi — 0. The other half 
of the real axis is the basin of attraction of a; = — oo. 



Phase portraits of the symplectic map (26 1 are de- 
picted in Figs. |3] and [4] and discussed in Subsection IIIC 
The singular line is located at a; = 1 and the focal point 
at z = (1,0). According to Eq. (37 1, the prefocal line is 
X — p + fi/e. The inverse (28 1 is not defined below this 
prefocal line for p > x — /i/e. 

An analytical calculation of the activation barrier can 
be carried out around the transcritical bifurcation at fi = 
1 by using Eq. ( 53 1 with the potential ( 50 1 yielding 



Wj? = |2^ - 2 - 2 In /i - (In ^)^ 



(62) 



which applies below and above the bifurcation. As for 
the logistic map, the activation barrier vanishes as the 
cube of the control parameter A/i = /.t — 1. This result is 
verified in Fig. [9] by comparison with the activation bar- 
rier computed by Monte-Carlo simulations. A very good 
agreement is also observed with the action (54 1 of the 
heteroclinic orbit h = W^P) n >V"(r). It show evidence 
that the heteroclinic orbits play the key role in the deter- 
mination of the activation barrier for the noise-induced 
escape from the attractor. 



FIG. 9: Action of the heteroclinic orbit Wh and activation 
barrier Wo (with error bars) from the Monte-Carlo simulation 
versus the parameter value fi for the exponential map ( 59 1 . 



The blue hatched curve is the analytical action calculated 
near the transcritical bifurcation and given by Eq. (62 1. 



B. Tangent bifurcation 

At the tangent bifurcation, a pair of stable and unsta- 
ble fixed points emerge in the system. A one-dimensional 
map illustrating this bifurcation is the quadratic map 



which has the fixed points 

Xi^2 = ±\/a* 



(63) 



(64) 



for II > 0. xi is stable while X2 is unstable. For > 0, 



the phase portrait of the symplectic map ( 26 1 is similar 



to the one of the logistic map but the elliptic island no 
longer exists for /i < where all the trajectories of the 
noiseless map quadmap escape to infinity. For /i > 0, the 
escape form the basin of attraction of Xi ^ x-p = 
should be activated by the noise. Close to the tangent 
bifurcation at = 0, the critical slowing down justifies 
the treatment with the continuous-time approximation. 
Using Eq. ( 53 1 with the potential ( 50 ) , the activation 



barrier is here given by 



3/2 



(65) 



Here, we recover the result obtained by Beale that the 
exponent 3/2 is universal for the activation barrier close 
to a tangent bifurcation [21 . 



into the logistic map ( 55 ) by the changes x 
V/4 



We notice that the quadratic map ( 63 ) is transformed 
/i)/2 and /i 



{fj, — i-Yl^- This transformation is fold- 
ing the parameter space two to one, therefore, mapping 
the transcritical bifurcation onto the tangent bifurcation. 
In this way, the exponent 3 of the transcritical bifurca- 
tion is consistent with the exponent 3/2 of the tangent 
bifurcation 121]. 
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C. Pitchfork bifurcation 

In the pitchfork bifurcation, a fixed point which pre- 
exists to the bifurcation destabihzes and generated two 
new stable fixed points. These two new attractors are 
separated by the now unstable fixed point, which is the 
boundary between the two basins of attraction. At the 
bifurcation, the two attractors meet with this boundary 
so that the pitchfork bifurcation is a boundary bifurca- 
tion in this regard. 

We are concerned by the first-exit process induced by 
the noise from one basin of attraction to the other across 
their common boundary, which constitutes a activation 
barrier between two wells for the kinetic potential ( 50 1 . 



As explained in Subsection |IIC| this first-exit process is 
similar to an escape process in the sense that its rate is of 
Arrhenius type and also given by Eq. (15) in the weak- 



noise limit. At the pitchfork bifurcation, we expect that 
the activation barrier Wq should vanish in a way charac- 
teristic of the bifurcation. Two maps are investigated to 
determine this barrier. 
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FIG. 10: Action of the heteroclinic orbit Wh and activation 
barrier Wo (with error bars) from the Monte-Carlo simulation 
versus the parameter value fi for the cubic map ( |66[ ). The blue 
hatched curve is the analytical action calculated near the 
pitchfork bifurcation and given by Eq. ( 69 1 . 



The cubic map 



Cubic map 



(66) 



undergoes a pitchfork bifurcation at the critical parame- 
ter value n = 1 and presents bistability above this bifur- 
cation for fj, > 1. For our purposes, we consider positive 
values of the parameter fi > 0. The cubic map has one 
fixed point below the bifurcation and three above: 



xi ^ 



if /i > 1 



(67) 
(68) 



Moreover, the cubic map has an unstable period-two or- 
bit at X4^5 = ±-y/l -t- 1/fi. For ii < 1, the fixed point 
xi = is an attractor and its basin of attraction is 
bordered by the period-two orbit: B{xi) =]x5,X4[. For 
/I > 1, the fixed point xi = becomes unstable and 
the fixed points X2 and xt, are two attractors with their 
basin of attraction given by B{x2) —]xi — 0,a;4[ and 
B{x^) —]xz,xi = 0[, respectively. Outside the inter- 
val delimited by the period-two orbit, trajectories can 
escape to infinity, which is a further attractor of the 
cubic map. The symplectic map (26 1 associated with 



the cubic map has three hyperbolic fixed points in the 
invariant subspace p = for /i > 1. Now, there ex- 
ist two elliptic islands symmetrically located around the 
origin. The origin /3 — (0, 0) stands for the barrier 
while r — (a;2,0) for instance. A heteroclinic orbit 
h = yV^(/?) n yV"(r) connects the hyperbolic fixed point 
r to /3. The associated symplectic map ( 26 ) has two focal 
points 21-1- = (±l/\/3, 0) and two corresponding prefocal 



For /i > 1, noise induces random jumps between the 
two basins Bix'i) and B{x-i), as well as escape toward 
infinity. The barrier between the two basins B{x-2) and 
B(x^ is located at the origin x\ = and vanishes at 
the pitchfork bifurcation while the barrier for escaping 
toward infinity is the period-two orbit, which does not 
vanish around [i — \. Accordingly, the escape events 
remain rare and negligible with respect to the first-exit 
events from one basin to the other. The difference of time 
scale between first-exit and escape events allows us to in- 
vestigate the former without perturbation by the latter 
in Monte-Carlo simulations in the parametric domain of 
interest. In the weak-noise limit, the first-exit process is 
controlled by the heteroclinic orbit h = n yV"(r). 

Using the continuous-time approximation, we can calcu- 
late analytically the activation barrier near the pitchfork 
bifurcation with Eq. ( 53 1 for the cubic map: 



(69) 



which holds for // > 1. The excellent agreement with the 



lines (37 1 



results of Monte-Carlo simulations is shown in Fig. 10 
where the action Wh of the heteroclinic orbit is also plot- 
ted. Here, the activation barrier vanishes as the square 
of the control parameter A/i = /i — 1. 



2. Gaussian map 

Another example of bistable system is given by the 
Gaussian map: 

a;„+i = /X exp(— a;^) with /i > (70) 
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FIG. 11: Action of the heteroclinic orbit Wh and activation 
barrier (with error bars) from the Monte-Carlo simulation 
versus the parameter value pu for the Gaussian map ( |70[ ). The 
blue hatched curve is the analytical action calculated near 
the pitchfork bifurcation and given by Eq. ( 73 1 . 



Its fixed points are 




(71) 
(72) 



The pitchfork bifurcation happens at /i = 1 and bista- 
bility manifests itself for /i > 1. Contrary to the cubic 
map, escape to infinity is not possible for this map, even 
in the presence of noise. In this regard, this model is very 
convenient to investigate the effects of noise on bistabil- 
ity. For < /i < 1, the whole real line is the basin of 
attraction of the origin x\ = 0. For fj, > 1, the basins of 
attractions of the two attractors emerging at the pitch- 
fork bifurcation are respectively B{x2) =] — oo,0[ and 
B{xz) =]0,-|-cx)[. Here also, the symplectic map (26 1 has 
two focal points z± ~ (±1 /\/2, 0) and two corresponding 



prefocal lines (37 1 



Taking xp — Q and xy = X2 for instance and calcu- 
lating the action of the heteroclinic orbit h — yV^(/3) n 
yV"(r) in the continuous-time approximation, we find 



/i — 1 — In yU 



(73) 



which holds for ^ > 1. Here again, the comparison with 
the results of Monte-Carlo simulations and the action 
of the symplectic map itself are excellent as shown in 
Fig. |11[ As for the cubic map, the activation barrier 
vanishes as the square of the control parameter A/z = 
/i — 1, which characterizes the pitchfork bifurcation. 



VI. CONCLUSIONS 

In this paper, we have studied noise-induced escape 
or exit from bifurcating fixed points of one-dimensional 
maps. Such fixed points correspond to periodic orbits 



of strongly dissipative dynamical systems. In this study, 
we have considered the transcritical, tangent, and pitch- 
fork bifurcations [37] . We followed a symplectic approach 
which applies in the weak-noise limit. For noise generated 
by independent Gaussian random variables, the time evo- 
lution can be expressed in terms of path integrals defined 
by an action functional. In the weak-noise limit, path 
integrals are dominated by the contributions of trajecto- 
ries which are the extremals of the action functional. For 
the present discrete-time systems, these trajectories are 
ruled by a symplectic map. This two-dimensional area- 
preserving map reduces to the noiseless one-dimensional 
map on its one-dimensional phase space which is left in- 
variant. 

The noninvertibility of the one-dimensional map and 
the associated symplectic map is at the origin of spe- 
cial phase-space structures. In particular, the two- 
dimensional map is not defined on singular lines, ex- 
cept at their intersections with the invariant subspace. 
These intersections are the so-called focal points |^ I43j . 
The multiple inverses of the symplectic map are de- 
fined on domains bordered by so-called prefocal lines 
33] • Strangely enough, the noninvertibility has for 
consequence that each focal point is the pre-image of a 
corresponding prefocal line. Under such circumstances, 
invariant curves such as the stable and unstable sets may 
form loops attached to the focal points and they may 
be composed of several disjoint pieces, as here observed 
for the global stable set. Otherwise, the symplectic map 
forms typical phase-space structures such as elliptic is- 
lands surrounded by chaotic zones and homoclinic or het- 
eroclinic tangles. Heteroclinic orbits have been identified 
which connect the hyperbolic fixed points corresponding 
to the attractor and the top of the barrier separating the 
basins of attraction. 

Noise induces the escape or exit from each basin of 
attraction, changing the attractors of the noiseless map 
into metastable states. The rates of escape or first exit 
have Arrhenius' dependence on the noise amplitude. In 
the present non-thermal stochastic systems, the analogue 
of the activation energy is given by the action of the het- 
eroclinic orbit connecting the attractor to the top of the 
barrier. We have developed an analytical method to cal- 
culate these actions close to the bifurcations. Thanks 
to the critical slowing down near the bifurcation, the 
noisy map behaves as a continuous-time stochastic sys- 
tem, whereupon the corresponding symplectic map can 
be approximated by a Hamiltonian flow with one degree 
of freedom [571 HH]- For this flow, the action of the hete- 
roclinic orbit can be calculated analytically. 

Applying this method to several one-dimensional noisy 
maps, we have shown how the activation barrier for 
noise-induced escape vanishes with the control parameter 
A/i = /i — /ic near the transcritical, tangent, and pitch- 
fork bifurcations. At these bifurcations, the activation 
barrier scales as Wq ~ |A/i|" with a universal exponent 
a characteristic of the bifurcation, as summarized in Ta- 
ble [l] The analytical results are in excellent agreement 
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bifurcation action exponent 

transcritical ~ I^mI'^ 3 

tangent ~ A/i'^/^ 3/2 

pitchfork ~ A/j^ 2 

TABLE I: Scaling behavior of the action Wo of the activation 
barrier versus the control parameter A/j, = jj, — fj,c and uni- 
versal exponent for noise-induced escape or first exit from the 

attractor undergoing different bifurcations. The correspond- 
ing rate is given by 7,7 ~ exp(— Wo/e) in terms of the noise 
amplitude e. 



with Monte-Carlo simulations. 

In conclusion, the symplectic approach is a powerful 



method to deal with stochastic dissipative systems in the 
weak-noise limit and to delineate the universality of their 
properties near bifurcations. 
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